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TECHNICAL PAPER 


ALGORITHM FOR ASTRONOMICAL, POINT SOURCE, SIGNAL 
TO NOISE RATIO CALCULATIONS 

I. GENERAL DESCRIPTION 


The starting point for this program is three other programs [1] (OTF, PSF and 
Camera) developed by Dr. D. J. Schroeder. References 2 through 4 provide the 
basis for the development work. The programs OTF and PSF compute the mono- 
chromatic point spread function for a star of a given visual magnitude as seen through 
a telescope with a circular obscuration ratio of 0.33, entrance pupil diameter D and 
system focal ratio F outside the Earth's atmosphere. In the absence of any degrada- 
tions, the point spread function (PSF) is computed from the modulation transfer 
function (MTF) of a perfect lens with a central obscuration. Image degradations can 
be introduced by multiplying the perfect lens MTF by the MTF's associated with the 
different types of degradations. The above programs include a MTF for image jitter, 
and MTFs for high, mid and low frequency wave aberrations associated with the 
optical telescope assembly. The PSF is assumed to be radiaUy symmetric to con- 
siderably simplify to the computations. No ray aberrations due to optical path dif- 
ferences, especially those associated with the off axis imaging of the telescope, are 
considered. Any off axis defects are assumed to be compensated for by the following 
instrument. The change in the PSF for a selected wavelength can be examined for 
different combinations and magnitudes of the degradations included in the program. 

Program Camera, in conjunction with OTF and PSF, adds a detector pel array, 
to the focal plane, a uniform cosmic background and the option to add a second star 
in the vicinity of the target star. For a given wavelength and for each pel in the 
detector array (maximum array size is 20 x 8) two relative signals are computed. 

The first is a target star relative signal. The second is the target star plus cosmic 
background signal. If a second star has been included, three relative signals are 
computed: the target star relative signal, the second star relative signal, and the 

target star plus second star plus cosmic background relative signal. 

The above program has been modified to also compute the encircled energy for 
the monochromatic PSFs, the polychromatic PSF, and signal to noise ratios (S/N) for 
a given wavelength interval as a function of observation time and observation times 
as a function of signal to noise ratios using the combined Space Telescope and Wide 
Field /Planetary Camera quantum efficiencies [5], filter functions [6], and a detector 
characterized by an rms readout noise, a mean dark current, and an optional dead 
space between the detector pels. 


II. USER INSTRUCTIONS 


The user is first asked to input the visual magnitude of the signal star and 
then to choose whether or not the spectral distribution of the star is uniform or 
blackbody. The cosmic background and background star (if requested) will eilso have 
the same type of distribution. If a blackbody distribution is chosen, the user is 
asked to input the effective temperature (degrees Kelvin) of the signal star. 



The effective temperatures of the cosmic background and background star (if 
requested) will be asked for later. 

Next the user is asked to input the number of wavelengths to consider and the 
lower and upper wavelengths (cm). For Space Telescope, the lower and upper wave- 
lengths are 10 ^ and 11 x 10 ^ cm. The maximum number of wavelengths that the 

program can handle is 21, which makes the wavelength increment 0.5 x 10~^ cm. 

For this case, the filter functions (0 < transmission < 1) and quantum efficiencies 
(electrons /photon) must be specified and read in (all 21 of them) in increments of 
“5 

0. 5 X 10 cm. The user is asked if there is a table of quantum efficiencies and 
filter functions. If the answer is yes to both questions then two files of data must 
have been previously prepared, which are now read into the program. One file con- 
tains the quantum efficiencies and the other contains the filter transmission. Because 
the program integrates over wavelength, the first and last non- zero filter trans- 
missions should be divided by 2 in order to use the trapezoidal rule for numerical 
integration. If the answer to both questions is no, the user is assumed to desire a 
system with 100 percent quantum efficiency and 100 percent transmission. 

To obtain a monochromatic star and background, the lower wavelength will be 
chosen as the desired wavelength, if the number of wavelength calculations is set to 

1. The user is then asked to input the quantum efficiency and filter transmission 
for that wavelength. 

Next the user is asked to input the diameter (cm) of the telescope entrance 
pupil, the focal ratio of the optical system and the central obscuration ratio (ratio 
of obscured to clear entrance pupil diameter) . For Space Telescope the obscuration 
Vatio is 0.33. 

The user is then asked to input the image degradation factors of jitter and 
high, mid and low frequency wave aberrations. For Space Telescope [1] the rms 
image jitter is anticipated to be 0.007 arc sec or less and the high, mid and low 

-5 ~5 

frequency aberrations are reported to be 0.121 x 10 cm, 0.1304 x 10 cm, and 
-5 

0.2361 X 10 cm, respectively. 

At this stage of the program the user may request a printout of the mono- 
chromatic point spread function and encircled energy. If this option is requested, 
the program outputs the Airy radius in arc seconds and cm using the shortest wave- 
length for which the product of the quantum efficiency and filter function is greater 
-4 

than 10 . The user is then asked to input the maximum radius (arc seconds) to be 

considered and the number of radius calculations desired. If the maximum radius to 
be considered is 0.2 and the number of radius calculations is 51, the point spread 
function and encircled energy values will be output in the range from 0 to 0.2 arc 
sec in increments of 0.004 arc sec. The wavelength is selected by entering a wave- 

~5 

length number from 1 to 21 for space telescope. For example, 9 is 5 x 10 cm, 10 
-5 -5 

is 5.5 X 10 cm, and 11 is 6 x 10 cm. The program then outputs the normalized 
point spread function [PSF(0)=1] , the encircled energy, EE, [EE(«>)=1], the nor- 

_2 -1 

malized point spread function scaled in units of energy cm sec as a function of 
radius along with the radius in units of cm and arc seconds. After the printout has 
been completed the user is asked if there is a desire to change any of the PSF 
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parameters. A yes answer allows the user to re-input a maximum radius, change the 
number of radius calculations or select another wavelength. The procedure described 
in this paragraph is repeated until a no answer is received for changing the PSF 
parameters. 

Next, the user is asked to enter the x,y widths of a pel (arc seconds) and the 
x,y pel center separation (arc seconds), the number of pels in the x,y directions and 
the x,y coordinates of the signal star chief ray in the detector plane. To visualize 
the effects of the inputs, see Figure 1, using the definitions listed below: 

XPW = X direction width of the CCD pel in arc seconds 

YPW = y direction width of the CCD pel in arc seconds 

XPD = X distance between pel centers in arc seconds, XPD > XPW 

YPD = y distance between pel centers in arc seconds, YPD > YPW 

NPX = number of pels in the x direction 

PSX = number of pels in the y direction 

PSX = X coordinate (arc seconds) of signal star chief ray, 0 < PSX < NPX. PD /2 

PSY = y coordinate (arc seconds) of signal star chief ray, 0 < PSY < NPY.YPD/2 

Dead space between the pels is created whenever XPD > XPW and/or YPD > YPW, 
i.e., the pel center separation is greater than the pel width in the same direction. 

The chief ray coordinates, which are now considered to be the optic axis, should be 
as close or closer to the origin than to the farthest corner of the pel that is farthest 
from origin. Otherwise, redundant situations are created that serve no purpose 
other than to complicate the programming. Smaller array sizes produce more accurate 
results in the (S/N) calculations because the PSF is computed at 200 different radii 
that range from zero to a value in arc seconds that is equal to the distance from the 
signal star chief ray coordinates to the coordinates of the farthest corner of the 
farthest pel from the origin (i.e., the larger the array size, the larger radius incre- 
ment) . The normalized point spread function is integrated over a pel to determine 
the fraction of energy, which will be denoted as diffraction efficiency, that falls 
within the' pel. This integration requires radial interpolation of the point spread 
function, which produces less error when the radial increments are made smaller. 

From these inputs just discussed in the paragraph, the program computes and outputs 
the diffraction efficiency as a function of wavelength for each pel in the array. 

Next, the user is asked to enter the visual magnitude of the cosmic background, 
and if a blackbody distribution was previously selected, an effective background tem- 
perature (degrees Kelvin) will be requested. The user is also asked whether or not 
to include a background star. If the answer is yes the user is asked for the visual 
magnitude of the background star, and if a blackbody distribution was previously 
selected, the effective temperature (degrees Kelvin) of the star will be requested. 

The x,y coordinates (arc seconds) of the background star chief ray are now 
requested. The diffraction efficiency, as discussed before, is computed for this 
star as a function of wavelength for each pel. In this case the point spread is also 
computed at 200 different radii. If the chief ray coordinates lie within the detector 
array, the radius ranges from a value of zero to a value equal to the distance from 
the chief ray coordinates to the farthest corner of the pel that is farthest from the 
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chief ray. If the chief ray coordinates lie outside the detector array, the radius 
ranges from a value that is equal to the minimum distance from the chief ray to the 
nearest pel to a value that is equal to the distance from the chief ray to the farthest 
corner of the farthest pel. From these just discussed inputs, the program computes 
and outputs the diffraction efficiency of the background star as a function of wave- 
length for each pel in the array. 

To complete the detector characteristics, the user inputs the rms readout noise 
in electrons per pel and the detector temperature in degrees Kelvin. For the par- 
ticvilar CCD array used with the Space Telescope and Wide Field /Planetary Camera, 
the readout noise has been suggested [5] to range from 13.9 to 17.8 electrons per 
pel rms. For the above CCD array, the maximum operating temperature [5] is 
expected to be around 178°K (-95°C), and the mean dark current in electrons per pel 
per second is computed from this output temperature. The user can input a mean 
dark current directly by first entering a temperature less than 4°K. A prompt will 
then appear requesting the desired dark current. 

The next set of inputs are concerned with the output data that is desired. 

The user first specifies the observation start and end times in seconds and the 
number of time calculations. For example, to compute the (S/N) in hour intervals 
from 1 to 14 hours, the user woiild input a start time of 3600 seconds, and an end 
time of 50400 seconds and 14 is the number of time calcxilations . To reverse the 
situation and determine the observation times needed to achieve a desired range of 
(S/N), the user inputs a start and end (S/N) and the desired number of (S/N) cal- 
culations. The program will proceed to output the desired data. 

After the output is complete, the user may repeat the calculations for a signal 
star with a different visual magnitude and/or a cosmic background with a different 
visual magnitude and/or a background star with a different visual magnitude, without 
having to re-input all the previous input data. An additional repeat calculation that 
can be accomplished with or without the magnitude changes, is to change the charac- 
teristics of the detector. This option, however, requires that the user re-input all 
data starting with the detector readout noise and ending with the number of (S/N) 
calculations. The program terminates when the user answers no to all the repeat 
calculation options. 


II. EQUATIONS USED 


The system of units used is the cgs system, except for the image jitter and 
detector focal plane coordinates and dimensions, which are in arc seconds. The total 
energy emitted in the eye responsive spectral region from a star of visual magnitude 
m and received per square meter per second outside the Earth's atmosphere is given 
by [7] 


Ij^ = 2.54 X 10 ® X 10“®-'*"’ lux 


( 1 ) 


The above equation is converted from photometric to radiometric units by dividing by 
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0.68 


9 


( 2 ) 


7.6x10 ^ cm 


/ 


-2 -1 -1 
K(X) dA lux per (ergs cm sec A A ) 


3.8x10"^ cm 


where A A is the wavelength interval (expressed in cm) and K(A) is the photopic eye 
response [7] given in Table 1. Dividing equation (1) by equation (2) is tantamount 
to assuming a uniform spectral distribution for the star and gives 


m 


3.7353xl0“®xl0 -2 ^ -1 . ,-l 

ergs cm sec A X 


I 


7.6x10 ^cm 
K(A) dA 
3.8x10 ^cm 


(3) 


If B(A,T) is the blackbody distribution having an effective temperature T (degrees 
Kelvin), then equation (3) can be converted to a blackbody distribution in the follow- 
ing manner [ 8] ; 


T rp^_ 3.7353x10 ®xl0 B(A,T) 

%(A,1) -g 

7.6x10 cm 


-2 -1 -1 
ergs cm sec A A 


/ 


K(A)B(A,T)dA 


3.8x10 ^cm 


(4) 


Any type distribution can be used for B(A,T) in equation (4), since multiplying equa- 

-5 

tion (4) by 0.68K(A) and integrating from 3.8 x 10 cm to 7.6 x 10 cm reproduces 
equation (1). If B(A,T) were a uniform distribution, then equation (3) would be 
obtained again. Assume that a uniform distribution has been selected and that the 
following quantities have been specified: 


2 

A = area of the telescope entrance pupil (cm ) 
he /A = ergs /photon 

Q(A) = quantum efficiency in electrons per photon 
F(A) = optical filter transmission 
DE (A) = diffraction efficiency 

The quantum efficiencies [5] used in the output example are shown in Table 2. This 
includes the combined effects of the optical telescope assembly and the Wide Field/ 
Planetary Camera. The filter function [6] used in the output example is listed in 
Table 3. The diffraction efficiency is qvute involved and will be discussed in 
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detail later. However, it is the integral of the normalized point spread function over 
the area of a pel. For a given wavelength, this integral gives the fraction of energy 
intercepted by that pel, since the integral over all space of the normalized point 
spread function is one. 

The rate at which electrons are liberated from a particular pel per wavelength 
interval at a chosen wavelength is given by 


®s^^^ ~ ^Ijjj(^)Q(^)F(X)DEg(X) electrons sec ^ AX 


(5) 


Equation (5) is called the monochromatic signal current. Integrating equation (5) 
over the wavelength region influenced by the quantum efficiency and optical filter 
produces the polychromatic signal current. If there is a dead space between the 
detector pels, the diffraction efficiency is adjusted to account for the dead space and 
the location at the chief ray of the diffraction pattern relative to the pel in question. 
The polychromatic signal current is given by 

A r ^ -1 

®s “ he J ^Ijjj(X)Q(X)F(X)DEg(X)dX electrons sec . (6) 


By the same analogy the signal current from a second star would be given by 
^2 

f electrons sec“^ , (7) 

where m' could be a different visual magnitude and DE^(X) is the diffraction efficiency 

of the second star for the same pel. The signal current for the cosmic background is 
given by 

^2 

^b " "^hc^^ f Inj’'0^)Q(^)F(^) electrons sec~^ , (8) 

where dQ is the area of a pel in arc seconds squared ( 113 ) . The cosmic background 
is treated as an extended source of visual magnitude m" and has units of 

per El , The dfi for an extended source replaces DE(X) for a point source. 

To compute a signal to noise ratio, the contributions to the signal term and the 
noise term must be defined. For the polychromatic case, for example, the signal is 
equation (6) multiplied by the time in seconds. If the incoming photon flux is 
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assumed to be Poisson distributed, then the variance of the signal is also equal to 
the mean signal. If the signal to noise ratio is defined as the ratio of the mean to 
the standard deviation of the mean, then the signal star will contribute to the noise. 
Other contributions to the noise include the cosmic background signal, the unwanted 
star signal if requested, the rms readout noise in electrons per pel and the mean 
dark ctirrent in electrons per pel per second. If all of the noise contributors are 
assumed to be statistically independent, then the signal to noise ratio is given by 


(S/N) = 




(S .t + S .t + S, .t + + D.t)^^^ 

sub 


(9) 


where R is the rms readout noise, D is the mean dark current and t is the time in 
seconds. If the second star option was not requested, would be set to zero. 

Equation (9) is computed for every pel and the only quantity that changes in the 

diffraction efficiency in the expressions for S and S . As previously mentioned, 

s u. 

the user has the option of inputing the dark current directly or inputing the detec- 
tor temperature and having the dark current computed for the CCD array on Space 
Telescope for the Wide Field /Planetary Camera. The equations (9) for the dark 
current are given by 


D = 29.3 X 10^ X T^'^^ exp(-5802.1Z) 


( 10 ) 


where 


1.1557 7.021 X 10 ^ T 

T 1108 + T 


( 11 ) 


and where T is the detector temperature in degrees Kelvin. 

Equation (9) can also be solved for the time required to achieve a particidar 
(S/N). Solving for the time gives 


t = 


(S/N)' 


(S +S +S^.+D) 
^ s u b ’ 


2S 


(S/N)^ 2 2 

^ ^ ^ + (S/N)^ (R/S^)^ 

4S ^ ® 


1/2 

sec . 
( 12 ) 


How the diffraction efficiency is calculated will now be discussed. As previously 
mentioned, the diffraction efficiency is the integral of the normalized point spread 
function over the detector pel area. The analytic expression for the normalized 
rotationally symmetric point spread function is given by [2] 
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(13) 


1 

PSF(r,X) = ^ f MTF(v^ 

P 0 


X) JQ(2irr [Dp. tt/X. 64800] v^) v^dv^ 


9 


where r is the radius of the PSF in arc seconds, D is the diameter of the telescope 

entrance pupil (cm), X is the wavelength (cm), is the normalized spatial frequency, 

is a zero order Bessel function of the first kind, ip is the pupil transmittance, 

648000 is the number of arc seconds in pi radians, and MTF is the modulation trans- 

fer function. For a central obstruction, the pupil transmittance is given by 1-e , 

where e is the ratio of the obstruction diameter to the clear diameter. Equation (13) 

is normalized such that PSF(0,A) = 1, however it can be scaled [2] in units of energy 

— 2—1 22 22 
cm sec by multiplying by t ir D /16A F , where F is the system focal ratio. 

Jr 

The MTF is composed of the product of the individual MTF's representing the perfect 
lens, the image jitter, and the high, mid and low frequency aberrations of the 
optical telescope assembly. The MTF equation for the perfect lens is given by [2] 

TPL(v^) = (A+B+O/d-e^) , (14) 

where 

A = I [arc cos(v^)-v^(l-Vj^2)^/^J 
A = 0 ; ^ 


B = ^ {arc cos(v„/e)-(v„/e)[l-(v„/e)^]^^^} ; 0 < v„/e < 1.0 

(16) 

B=0 ; v^/e>1.0 


0 ^ "n < 


1.0 


(15) 


C = -2e^ ; 0 1 < (l-e)/2 

C = ^ |e sin<j) + ^ (1+e^) - (1-e^) arc tan tan j ~ 2e^ 

(l-e/2) < < (l+s/2) 


C = 0 ; > (l+e/2) 


(17) 


and 


(j) = arc cos 


’i 2 . 2 

l-e -4v^ 
2e 
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The MTF for image jitter is given by [2] 


TJ(v^,A) = exp 


„ 2 . 2 . 2 . 2 . 2 
^3 % ^ % 

• 648000^ 


( 18 ) 


where is the rms jitter in arc seconds, Dp is the telescope entrance pupil diameter 

(cm) , 648000 is the number of arc seconds in pi radians and A is the wavelength 
(cm). The MTF for the high frequency wave aberrations is given by [4] 


TH(A) = exp 



(19) 


where is the rms high frequency aberration (cm). The MTF for the mid frequency 
aberrations are given by [1] 



< 1/18 (20) 


where is the rms mid frequency aberration (cm) . The MTF for the low frequency 
aberrations is given by [ 1] 


2.815 o-r / 5irv \ 

TL(v^,A) = 1 ^ i^sin(-^l ; < 0.355 (21) 

Or 

TL(v^,A) 1 --oToayrr sec[3/4(Vj^-0.4)D ; > 0.355 

where cr2is the rms low frequency aberration (cm). If TL(v^,A) becomes negative, 
it is set to zero. Thus, the MTF(v^^,A) in equation (13) is given by 

MTF(v^,A) = TPL(v^) .TJ(v^,A) .TH(A) .TM(v^^,A) .TL(v^,A) (22) 

Fora given wavelength, the encircled energy, EE(r^,A), is given by [2] 


9 





EE(r^,X) 


648000 


2it r 


I f 


PSF(r,X,0) rdrde 


(23) 


Since PSF is radially symmetric, equation (23) becomes 


II^T 

EE(r^.X) = -2^ 


°p • ' 


X • 648000 


r ° 

I PSF(r,X)rdr 
o 


(24) 


The encircled energy is normalized such that EE(oo,X) = 1. Substituting (13) into 
(24) gives 


EE(r^,X) = 4n 


or 


EE(r^,X) = 2iTr^ 


Dp • . 

X • 648000 


-1 2 1 


J MTF(v^,X) 


/ ■ 648000 ) 

- o ^ 




( 25 ) 


Dp . . 

X • 648000 


/ r °p * ’ 

J MTF(Vp.x) Jj(2Wp 

n ^ 


648000 


V ) dv 
n'^ n 


where is a first order Bessel function of the first kind, is the diameter (cm) 

of the telescope entrance pupil, r^ the radius of the encircled energy in arc seconds, 

X is the wavelength and 648000 is the number of arc seconds in pi radians. To 
achieve reasonable accuracy needs to be calculated for at least 400 points in the 

range 0 < < 1. Thus, is incremented in increments of 1/400. It is also 

recommended [1] that the following numerical integration scheme [10[ be used: 


b 

/ = A 

a 


1 . 4y 6 . 4y j^+ 2 . 4y 2 + 6 . 4y 3 + 2 . 8y ^ 


+6 . 4y g+2 . 4y g+6 . 4y ,^+2 . 8y g 


(26) 


where 4nh=b-a. 


] 


4n- 3'^^ • 4y .^+6 . 4y ^ +1 . 4y 


4n-2 


4n-l 


4n 
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The computation of the diffraction efficiency is very similar to computing 
encircled energy except that the integration is performed over a square or rectangu- 
lar area instead of a circular area. Thus, for a particiilar pel, equation (23) becomes 


ITT. 


DE(A) = 


r r 

I A • 648000 J J 


PSF(x,y,A) dxdy 


(27) 


where x,y are expressed in arc seconds. Each pel is broken into a 20 by 20 sub- 
pixel array. The distance from each subpixel to the chief ray of the diffraction 
image is computed. Since the PSF is radially symmetric, the value of PSF(x,y,A) is 

2 2 1/2 

found by interpolating the value of PS (r,A) to a value of r = (x +y ) . This 

procedure is used for the signal star as well as the second star. 
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TABLE 1. PHOTOPIC EYE RESPONSE 


X (10 ^cm) 

K(X) 

X (10 ®cm) 

K(X) 

X (10 ^cm) 

K(X) 

3.8 

0.00004 

5.1 

0.503 

6.4 

0.175 

3.9 

0.00012 

5.2 

0.71 

6.5 

0.107 

4.0 

0.0004 

5.3 

0.862 

6.6 

0.061 

4.1 

0.0012 

5.4 

0.954 

6.7 

0.032 

4.2 

0.004 

5.5 

0.995 

6.8 

0.017 

4.3 

0.0116 

5.6 

0.995 

6.9 

0.0082 

4.4 

0.023 

5.7 

0.952 

7.0 

0.0041 

4.5 

0.038 

5.8 

0.87 

7.1 

0.0021 

4.7 

0.09 

6.0 

0.631 

7.3 

0.00051 

4.8 

0.139 

6.1 

0.503 

7.4 

0.00025 

4.9 

0.208 

6.2 

0.381 

7.5 

0.00012 

5.0 

0.323 

6.3 

0.265 

7.6 

0.00006 



TABLE 

2. QUANTUM EFFICIENCIES 


o 

3 

Q(X) 

X (10 ®cm) 

Q(X) 

X (10 ®cm) 

Q(X) 

1.0 

0.0^ 

4.5 

0.08^ 

8.0 

0.12 

1.5 

0.028^ 

5.0 

0.12 

8.5 

0.105^ 

2.0 

0.05 

5.5 

0.15^ 

9.0 

0.09 

2.5 

0.05 

6.0 

0.18 

9.5 

0.08^ 

3.0 

0.05 

6.5 

0.18^ 

10.0 

0.07 

3.5 

0.06 

7.0 

0.18 

10.5 

0.045* 

4.0 

0.04 

7.5 

0.15^ 

11.0 

0.02 


♦Interpolated Values 

TABLE 3. FILTER FUNCTION 


(10 ® cm) 

F(X) 

5.0 

0.278* 

5.5 

0.899 

6.0 

0.578 

6.5 

0.0835* 


♦The values were divided by 2 for using trapezoidal 
rule in numerical integration over wavelength. 
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APPENDIX A 


OUTPUT EXAMPLE 


ENTER SIGNAL STAR(tl) UlSUAL MAGHHUDE 
727 , 

ENTER CHOICE OF SPECTRAL DISTRIBUTION FOR SOURCE 

FLAT DI$TRIBUTION*e|BLACICBODy DISTRIBUTION-1 

? 0 , 

ENTER NUNBER OF UAUELCNGTHS 

? 21 , 

ENTER LOUER AND UPPER UAMELENGTHSCCN ) 

?KE-5ai.E-5 

IS THERE A TABLE OF OUANTUN EFFICIENCIES?< VES-l|NO-B ) 


IS THERE A TABLE OF FILTER TRANSmSSI0NS7(VES-liN0-6) 

? 1 , 

ENTER TELESCOPE ENTRANCE PUPIL DIANETERCCH) 

7240, 

ENTER SYSTEH FOCAL RATIO 
712.9 

ENTER OBSCURATION RATIO 
?e*33 

ENTER RHS IflAGE JITTER(ARC SECONDS) 

?e«ee7 

ENTER HIGH FREQUENCY RNS ABERATIONS(CA) 

79,121E-5 

ENTER niD FREQUENCY RNS ABERATI0N$(CN) 

7O.1304E-5 

ENTER LOU FREQUENCY RNS ABERAT10H$(CN> 

7e.2361C-5 

DO YOU WISH TO CONFUTE NORNALIZED POINT SPREAD 
FUNCTION,PSF(0)-1,AND THE NORNALIZED ENCIRCLED ENERGY, 
EE(INF)-l?(VES-liNO-0) 

71. 

NIHINUN AIRY RADIUS- .7869000^-*03 CN- .S242SS39E-01 ARC SECONDS 

ENTER NAXinUN RADIUS TO CONSIDERCARC SECONDS) 

70,5 

ENTER NUNBER OF RADIUS CALCULATIONS 
751. 

ENTER DESIRED UAUELENGTH NUNBERC 1,2, , . . ,21 ) 

79, 

CLEAR SCREEN AND RETURN TO CONTINUE 
7 
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PIAXIMUC' RADIUS* .S00«00ME 00 ARC SECONDS* . 7S0'49lS7E-08 CM 

HUmtP Of RADIUS CALCULATIONS St J UAUELENOTH* .S0000000E-04 CM 

RAMUS(Cn) PNT. SPREAD FN. ENCIRCLED ENERCV (ARC SECONDS) ENERCV/SEC/CH««8 
•0O00O000E 00 .71783407E 00 .00000000E 00 .00000OO0E 00 .5403O8S6E 11 

. 1S009831E-03 .62818690E 00 .80O34691E-01 . 100O0O0OE-01 .47807730E 11 

.30019663E-O3 .413S6896E 00 .8638S7S0C 00 .20000000E-01 .31474380E 11 

. 4S029494E-03 . 19860445E 00 .4375081 IE 00 .3000000OE-O1 .14668089E 11 

.60939335E-03 .57680580E-O1 .53230129E 00 . 40O0000OE-01 .43897407E 10 

.7S049157E-03 .20511074E-O1 .S6671843E 00 .50O00000E-O1 .15fi09811E 10 

.900S8988E-03 .32076381E-O1 .59954483E 00 .fi0O0O00OE-01 .24411S08E 10 

.1050S882E-O2 .41359947E-01 .65898910C 00 . 7000O00OE-01 .31476702E 10 

.12007865E-O2 .32244479E-01 .72702312E 00 .8OO0O0OOE-01 .24S39438E 10 

.13508848E-02 . 16218223E-01 .77586572E 00 . 9000OO00E-O1 .12341244E 10 

. 15O09831E-O2 .613438S4E-02 .79957848E 00 .1O0O00O0E 00 .46685316E 09 

.16510814E-02 .32642920E-02 .81042938E 00 .11OO00O0E 00 .24842668E 09 

.18911798E-02 .27686108E-02 .8184190SE 00 .1200OOOOE 00 .21070321E 09 

.19512781E-02 .21322428E-02 .82599344E 00 .13O0O00OE 00 .16227286E 09 

.21013764E-02 . 18672712E-02 .83234957E 00 .14OOO000E 00 .14210739E 09 

.22514747E-02 .211972O6E-02 .83921188E 00 .1SOO000OE 00 .16131987E 09 

.24015730E-02 . 196fi4822E-02 .84735678E 00 .1600000OE 00 .14965777E 09 

.25516713E-02 . 12263025E-O2 .85378890E 00 .170O00OOE 00 .93326904E 08 

.27017696E-02 . 83088452E-03 .85760O40E 00 .180000OOE 00 .83233892E 08 

.28518679E-02 . 12570514E-O2 .86245988E 00 .19O0O000E 00 .9S667029E 08 

.30919663E-02 . 1 794876SE-02 .86935760E 00 .290000O0E 00 .136S9784E 09 

.31520646E-O2 . 16565774E-02 .8781fiS18E 00 .21OO0OO0E 00 .12607268E 09 

.33021629E-02 . 1O036829E-O2 .88541796E 00 .22O000O0E 00 .76384S91E 08 

.34S22612E-02 .49217134E-03 .8894207SE OO .23OO0OO0E 00 .37458360E 08 

.36O23S95E-08 .34663041E-O3 .8917125SE 00 .240O00OOE 00 .2638O068E 08 

.37524578E-02 .326687672-03 .893SO920E 00 .25O00OOOE 00 .24862339E 08 

.39O25561E-02 .28592982E-03 .89511280E 00 .260000002 00 .217604912 08 

.40526545E-02 . 30576034E-03 .89710784E 00 .27000000E 00 .23269679E 08 

.42027528E-02 . 38363752E-03 .89876181E 00 .280000002 00 .29196469E 08 

.43528511E-02 .3712616SE-03 .90232322E 00 .29OOOO0OE 00 .28254611E 08 

.4S029494E-02 .25106O76E-O3 .90436156E 00 .30000000E 00 .19106805E 08 

.46530477E-02 .20103978E-03 .90672520E 00 .31000000E 00 .15299993E 08 

.48O31460E-O2 .312688942-03 .9081232SE 00 .32000000E 00 .23796975E 08 

.495324432-02 .437818172-03 .911959252 00 .330000002 00 .333198492 08 

.510334262-02 .408594522-03 .914753362 00 .340000002 00 .310958042 08 

.525344102-02 .263349912-03 .918096052 00 .350000002 00 .200420632 08 

.540353932-02 .152121662-03 .918847142 00 .360000002 00 .115771142 08 

.555363762-02 .122717942-03 .920503782 00 .370000002 00 .933936432 07 

.570373592-02 .116364232-03 .920551902 00 .380000002 00 .885581922 07 

.585383422-02 .105731962-03 .921735902 00 .390000002 00 .804589702 07 

.600393252-02 .118503552-03 .924090512 00 .400000002 00 .901863052 07 

.615403082-02 .146637072-03 .923949832 00 .410000002 00 .111597132 08 

.630412922-02 .140063522-03 .925567432 00 .420000002 00 .106594372 08 

.645422752-02 .991186862-04 .928019152 00 .430000002 00 .754335912 07 

.660432582-02 .879500772-04 .927317322 00 .440000002 00 .669337982 07 

.675442412-02 .133430532-03 .929505532 00 .450000002 00 .101546382 08 

.690452242-02 .178117962-03 .931031752 00 .460000002 00 .135555452 08 

.705462072-02 .163394232-03 .934239582 00 .470000002 00 .124350052 08 

.720471902-02 .108487082-03 .933661662 00 .480000002 00 .825633392 07 

.735481732-02 .689563112-04 .935599262 00 .490000002 00 .524787242 07 

.750491572-02 .587744632-04 .936120362 00 .500000002 00 .447298992 07 

DO YOU UISH TO CHANGE ANV PSF PARAn2TERS?(VES«ljNO*0) 

70. 

CLEAR SCREEN AND RETURN TO CONTINUE 
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ENTER X,Y PEL UIDTHSCRRC SECONDS) 

? 0 . 1 , 0.1 

ENTER X,V PEL CENTER $EPRRRTION$(RRC SECONDS) 

90 . 1 . 0 . 1 

ENTER NUNBER OF PELS IN X,V DIRECTIONS 

92 . , 2 . 


ENTER X,V COORDINATES OF SIGNAL STAR (ARC SECONDS) 


7e.es 

UAU- 

.^^eeeeeE-e4 

rx-i 

IV- 

1 

DIFE- 

.57745e45E $$ 

UAU- 

*seeeeeeeE-e4 

IX-l 

IV- 

8 

Dire- 

p4794ei38E*ei 

UAM- 

.seeeeeeeE-e4 

IX-8 

IV- 

1 

OIFE- 

•4794ei38E-ei 

UAM- 

.seeeeeeeE-#4 

ix-a 

IV- 

8 

DIFE- 

«ie863ei4E-ei 

UAU« 

.sseeeeeeE-e4 

IX-l 

IV- 

1 

OIFE- 

o57873858E 88 

UAU- 

*s5eeeeeeE-e4 

IX-l 

IV- 

8 

DIFE- 

p47984166£-81 

UAU- 

.sseeeeeeE-e4 

IX-8 

IV- 

1 

DIFE- 

.47984166E-81 

UAU* 

.sseeeeeee-e4 

IX-8 

IV- 

8 

DIFE- 

ol8838871E-81 

IMU- 

.6eeeeeeec-e4 

IX-l 

IV- 

1 

DIFE- 

a5€934979E 88 


.6eeeeeeeE-e4 

IX-l 

IV- 

8 

DIFE- 

a46483964E-81 

UAU- 

.6eeeeeeee-e4 

IX-8 

IV- 

1 

DIFE- 

a46483964E*81 

UAV- 

.6eeeeeee£-e4 

IX-8 

IV- 

8 

DXFE- 

*14751487E-81 

MAM- 

•6seeeeeeE-e4 

IX-l 

IV- 

1 

D2FE- 

*58484891E 88 

UAM- 

.65eeeeeeE-e4 

IX-l 

IV- 

8 

DXFE- 

•44S73175E-81 

MAU- 

.656eeeeeE-e4 

IX-8 

IV- 

1 

OIFE- 

o4457317SE-81 

UAU- 

.65eeeeee£-e4 

IX-8 

IV- 

8 

DIFE- 

a17411538£-81 


ENTER COSNIC BACKGROUND UISUAL MAGNITUDE 
923 . 


DO YOU WISH TO INCLUDE A BACKGROUND SrAR7(VES-t)N0-S> 


ENTER RNS READOUT HOISE(ELECTRONS/PEL) 

918 . 

ENTER DETECTOR TEHPERATURE(KELUIN) 

9178 . 

ENTER OBSERUATION START AND END TINES(SECONDS) 

93600 . . 50400 . 

ENTER NUNBEROl) OF TIME CALCULATIONS 
914 . 

ENTER START AND END SIGNAL TO NOISE RATIOS (S/N) 

91 . . 10 . 

ENTER NUNBEROl) OF (S/N) CAUUIATIONS 
910 . 

CLEAR SCREEN AND RETURN TO CONTINUE 
9 
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00 


TEieSCOPE ENTRANCE PUPIL DIANETER 849 .AAe Cfl 
OPTICAL SVSTEn FOCAL RATIO ia«99 

SOURCE STAR CHARACTERISTICS 
SOURCE NACNITUDE 87*99 

COSRIC BACKGROUND CHARACTERISTICS 

COSniC BACKGROUND MAGNITUDE 83*99 

SYSTEM SPECTRAL CHARACTERISTICS 

UAUELENGTH(CH) 0UAN*EFF1C* FILTER FUNCTION 
*5999E*94 *1899E 99 *8789E 99 

*5S99E-94 *1599E 99 *8999£ 99 

*6999E-94 *1899E 99 *5789E 99 

*6S99E-94 .1899E 69 *8359E*91 

DETECTCN^ CHARACTERISTICS 

PEL DINENSIONSCARC SECONDS )tX« *1699 V- *1969 

PEL CENTER SEPARAT10N(ARC SECONDS )IX» *1999 Y- 

SOURCE STAR COORDINATES (ARC SECONDS)tX« *9599 Y* 

ARRAY SIZEtX- 8 BY V* 8 

RNS READOUT N01S£(EL£CTR0HS/PEL >118*99 

DETECTOR TENPERATURE(KELUIN) 1178*9969 

DARK CURRENT ( ELECTRONS/PEL/SECOND ) i * 53387884E-98 

SIGNAL TO NOISE RATIO AS A FUNCTION Of OBSERUATION TINE 
AND PEL NUMBER 


TIME(SECONDS) 

(S/N) 

PEL 

NUMBER 

*36999E 

94 

•7986E 61 

1 

1 


94 

*18888 68 

1 

1 

*19896E 

95 

*16588 62 

1 

1 

*14496E 

65 

.19758 68 

1 

1 

*18969E 

95 

*88488 98 

1 

1 

*81696E 

95 

*84938 92 

1 

1 

*85896E 

95 

*87188 68 

1 

1 

*88896E 

95 

*29858 62 

1 

1 

*38496E 

95 

*31198 68 

1 

1 


95 

*33638 62 

1 

1 

*39696E 

95 

*3476E 68 

1 

1 

*43896E 

95 

*36488 62 

1 

1 

*46899E 

95 

*38668 62 

1 

1 

*59496£ 

95 

*39528 62 

1 

1 


94 

*77568 66 

1 

8 

*78696E 

94 

*13298 61 

1 

2 

*16896£ 

95 

.17788 61 

1 

8 

.14496E 

95 

*81478 61 

1 

2 

*18996£ 

95 

*84788 61 

1 

8 

*81696E 

95 

*87758 61 

1 

8 

*85296E 

95 

*36478 61 

1 

8 

*88896£ 

95 

*32998 61 

1 

8 

*384998 

95 

*35358 61 

1 

8 

*35996£ 

95 

*37578 61 

1 

8 

*39696E 

95 

*39688 91 

1 

8 

*438668 6S 

.41688 61 

1 

8 

«46866E 

95 

*43668 61 

1 

8 

*56496E 

65 

.4544E 61 

1 

8 

*359668 

94 

*77568 69 

8 

1 

*789668 

94 

*13898 61 

8 

1 

* 16896E 

95 

.17788 61 

8 

1 

*144968 

95 

.81478 61 

8 

1 

* 18996E 

65 

*84788 61 

2 

1 

*816968 

65 

*87758 61 

8 

1 


8S299E 

95 

•3947E 91 

2 

1 

888998 

95 

.32998 91 

8 

1 

32499E 

95 

.3535E 91 

2 

1 

36996E 

95 

.3757E 91 

8 

1 

396998 

95 

•3968E 91 

2 

1 

43869E 

95 

.4168E 91 

8 

1 

468998 

95 

.4369E 91 

8 

1 

59496E 

95 

.4544E 91 

8 

1 

369998 

94 

*2814E 96 

8 

2 

78999E 

94 

*38168 96 

8 

8 

198968 

95 

*51958 69 

2 

8 

144968 

95 

*6292E 96 

8 

8 

189668 

95 

.7168E 96 

8 

2 

.816698 

95 

.8938E 96 

8 

2 

858968 

65 

*88348 66 

8 

2 

2889^ 

95 

.9572E 96 

2 

2 

32496E 

95 

.19868 91 

8 

2 

369698 

65 

*19918 91 

2 

2 

396668 

95 

*11538 91 

2 

2 

432968 

95 

*12118 61 

8 

2 

468668 

95 

.1868E 91 

2 

2 

564968 

65 

.13818 61 

8 

2 


<S/H) 

SECONDS 


PEL 

NUMBER 

1.66 

*889498 

92 

1 

1 

2*96 

*817428 

63 

1 

1 

3.66 

.39891E 

63 

1 

1 

4*99 

*686688 

63 

1 

1 

5*96 

•99467E 

63 

1 

1 

6*96 

.184518 

64 

1 

1 

7*66 

.164498 

64 

1 

1 

8*96 

*819488 

64 

1 

1 

9*66 

*868348 

64 

1 

1 

16*66 

*389288 

94 

1 

1 

1*66 

.247698 

64 

1 

2 

2*66 

•99138E 

64 

1 

2 

3*66 

.198658 

65 

1 

2 

4*96 

.356498 

65 

1 

2 

5*66 

*545698 

65 

1 

2 

6*96 

.784268 

65 

1 

2 

7*66 

. 10668E 

96 

1 

2 

8*96 

*139158 

66 

1 

2 

9*66 

*176628 

66 

1 

2 

19*96 

.21783E 

66 

1 

2 

1*96 

*847698 

64 

2 

1 

2*66 

.991388 

64 

2 

1 

3*96 

*198658 

65 

8 

1 

4*96 

.356498 

65 

8 

1 

5*66 

.545698 

65 

2 

1 

6*96 

*784868 

65 

2 

1 

7*66 

*196628 

96 

2 

1 

8*96 

.139158 

96 

8 

1 

9*66 

•176628 

96 

2 

1 

19*96 

.217838 

96 

2 

1 

1*66 

*858948 

95 

2 

8 

8*96 

*198168 

96 

2 

2 

3*66 

.829868 

96 

2 

8 

4*96 

*497898 

96 

8 

2 

5*96 

.636188 

96 

8 

2 

6*96 

.91SS3E 96 

8 

2 


7. M 

8. eo 

9. ee 

ie«ee 


«ia465E •? 
aeaseE •? 
.aee^^e 8? 

.aS436E 87 


DO YOU UtSH TO CHANGE SOURCE NAGNITUDE?<VES-1|NO-0) 

70* 

DO YOU WISH TO CHANGE COSNIC BACKGROUND nAGNXTUDETCVES-iiNO-G) 

?e. 


DO YOU UISH TO CHANGE BACKGROUND STAR NAGNXTUD£7(YES*1)N0«G) 
70* 


DO YOU UISH TO CHANGE DETECTOR ErFECTS?(YES*l,N0*6) 
70* 

DID YOU MAKE ANY CHANGES?(VES«1|NO*0 ) 

70. 

TSTOPt 0 
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APPENDIX B 


PROGRAM LISTING 

07S^I6 AUG 10»'84 OC/PNTSRC.JAYROE 

IMPLICIT DOUBLE PRECISI0N(A-H,0-Z) 

DIMENSION TPL(400),TJ(400),THC2U 

DIMENSION TM(400»21),TL(a00,2n,FXl («00 ) , DFE ( 1 0, 1 0) 

DIMENSION SIGNL(200),PELSa0,10),PEL8(10/10),DIFE(10,10) 
COMMON /8ESF/FXO(«00)#PR(200,2Uf TTOT(«00,2n,OUEF(21)rFF(2l) 
COMMON /PHOF/XKC40) 

WRITEU02,6000) 

6000 FORMAK'ENTER SIGNAL STAR(#1) VISUAL MAGNITUDE') 
READ(l0l,i4000)XSl 

aOOO FORMATCaOG) 

XISla0.37353*10,0A*(-0,4*XSl-5.) 

HCal.R865E-l6 

HC0Kal,a388 

DLAal,E-6 

KRESPsO 

WRITEn02,6001) 

6001 FORMAT!/, 'ENTER CHOICE OF SPECTRAL DISTRIBUTION FOR SOURCE') 
WRITE(102,6002) 

6002 FORMAT('FLAT DISTRIBUTIONsO jBLACKBOOY DISTRIBUTIONS!') 
REAO(101,4000)KRESP 

IF(KRESP,E0.1)G0 TO 90 
XISLls(XISl*10,0**5)/(l,0682*HC) 

GO TO 120 

90 WRITE(102,6004) 

6004 FORMAT!/, 'ENTER SIGNAL STAR EFFECTIVE TEMPERATURE !KELVIN) ' ) 
REAO!101,4000)ST1 
DO no IKsl,39 
READ!20,5000)XK!IK) 

5000 FORMAT!E15,8) 
no CONTINUE 

CALL PH0I!ST1,HC0K,DLA,HC,XINT) 

XI8L1=XIS1*XINT 
120 WRITE!102,6008) 

6008 FORMAT!/, 'ENTER NUMBER OF WAVELENGTHS') 

READ!101,4000)NWAV 

DLsO, 

DO 125 lsl,21 
QUEFCDsl. 

FF!nsl, 

125 CONTINUE 
FF(l)s0.5 
FF!21)s0.5 

IF(NWAV.E0,1)60 TO 130 

WRITE(102,6012) 

6012 FORMAT!/, 'ENTER LOWER AND UPPER WAVELENGTHS !CM) ' ) 
REA0C101,4000)WAVL#WAVH 
XsNWAV-1 

DLs!WAVH«WAVL)/X 
GO TO 140 

130 WRITE(102^6016) 
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6016 FORMAK/, 'ENTER WAVELENGTH (CM) ' ) 

READ(101,^000)WAVL 

WRITE(102,6020) 

6020 FORMATC/, 'ENTER QUANTUM EFFICIENCY AND FILTER TRANSMISSION') 
REAOn01,4000)QUEF(l),FF(l) 

60 TO 180 

lao WRITE(102,602A) 

602« FORMAT(/,'IS THERE A TABLE OF QUANTUM EFF IC lENC lES? ( YES=1 ; N0=0 ) ' ) 
READ(101,4000)IRESP 
IF(IRESP,NE,1)G0 TO 160 
DO 150 I=1,NWAV 
150 REA0(22,5000)QUEF(I) 

160 WRITE(102,6028) 

6028 FORMAT!/, 'IS THERE A TABLE OF FILTER TRANSMISSI0NS?(YES=1 ;N0 = 0) ') 
REA0(101,«000)IRESP 
IFdRESP.NE.DGO TO 180 
DO 170 I=1,NWAV 
170 REA0(2«,5000)FF(I) 

180 WRITE(102,6032) 

6032 FORMAT!/, 'ENTER TELESCOPE ENTRANCE PUPIL 01 AMETER (CM) ' ) 
REA0!101,A000)DP 
WRITE(102,6036) 

6036 FORMAT!/, 'ENTER SYSTEM FOCAL RATIO') 

REA0(101,^O00)FR 

CUTF=1./!FR*WAVL) 

WRITE!102,6000) 

6000 FORMAT!/, 'ENTER OBSCURATION RATIO') 

READ(101,4000)EN 

EN2=EN*EN 

ENPLS=1,+EN2 

ENMNSsl,-EN2 

ENRAT=!1.+EN)/(1,-EN) 

PI=3, 1015926 
T0PI=2./PI 
DM=0, 0025063 
XARGs-DM 
00 210 1=1,400 
A = 0. 

8 = 0 , 

C=0, 

XARG=XARG+DM 
IF!XARG,GT.1,)G0 TO 210 

A»T0PI*(AC0S(XARG)-XARG*SQRT!1.-XARG*XARG)) 

IF!EN,LT. 0,001)60 TO 210 
YARG^XARG/EN 
IF!YARG,GT,1,)G0 TO 195 

B=T0PI*EN2*!AC0S(YARG)-YARG*3QRT!l,-YAR6*YARG)) 

195 2ARG=XARG*2./!1,-EN) 

IF(2ARG.GT,1,)G0 to 200 

C=-2,*EN2 

60 TO 210 

200 WARG«XARG*2,/!1,+EN) 
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IF(WARG,GT.1.)G0 TO 310 
PHI=AC0S(U.+EN2-^}.*XARG*XARG)/(2.*EN)) 
CsEN* 8IN(PHI)+PHI*ENPLS/2,-ENMNS*ATAN(ENRAT*TAN(PHI/2,)) 
C=C*T0PI-2,*EN2 
210 TPLn)s(A+B+C)/ENMNS 
ON=0,0035063*WAVL 
WRITE(102,604A) 

60AA FORMATC/, 'ENTER RMS IMAGE JITTERCARC SECONDS)') 
READ(101,A000)SIGJ 
WRITE(102,6Q«8) 

60A8 FORMATC/, 'ENTER HIGH FREQUENCY RMS ABERATIONS (CM) ' ) 
REA0U01,4000)SIGH 
^RITE(102,6052) 

6052 FORMAT(/#'ENTER MID FREQUENCY RMS ABERATIONS (CM) ' ) 
REA0(101,A000)SIGM 
WRITE(102,6056) 

6056 FORMAT!/, 'ENTER LOW FREQUENCY RMS ABERATIONS(CM) ' ) 
READ(101,^000)SIGU 
WAV=WAVL-DL 

FACJ9PI*SIGJ*DP*PI/6A8000, 

FACJs-2,*FACJ*FACJ 

FACHs-4.*PI*PI*SIGH^SIGH 

FACM=-A,*PI*PI*SIGM*SIGM 

DNU=DM 

DO 250 NW=1,NWAV 
PROD=OUEF(NW)*FF(NW) 

IA(AVsWAV + DL 
XNU=-DNU 

TH(NW)=EXP(FACH/WAV**2) 

DO 250 1=1,400 
XNUsXNU+DNU 

TJ(I)=EXP(FACJ*XNU*XNU/WAV**2) 

TT0T(I,NW)=0, 

IF(PR00.LT.l.E-4)60 TO 250 
XARGsXNU 

TM(I,NW)=EXP(FACM/WAV**2) 

YARG=l.-l8,<kXARG 
IFCYARG.lt. 0,)60 TO 220 
ZARG=-FACM*YARG**1.5/WAV**2 
TM(I,NW)=TM(I,NW)*EXP(ZARG) 

220 IF(XARG.GT.O,355)GO TO 230 
YARG*PI*XAR6/0.6 

TL(I,NW)s1,-2,815*SIGL*8IN(YARG)/WAV 
GO TO 240 

230 YARG=0.75*(XARG-0.4) 

TL(I,NW)sl.-SI6L*(-0.9+l./COS(YARG))/(0.0373*WAV) 

240 IF(TL(1,NW).LT.0.)TL(I,NW)=0. 

TT0T(I,NW)=TPL(I)*TJ(I)*TM(I,NW)*TL(I,NW)*THCNW) 

250 CONTINUE 

WRITEC102,6060) 

6060 FORMAT!/, 'DO YOU WISH TO COMPUTE NORMALIZED POINT SPREAD') 
WRITE(102,6064) 
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6064 FORMATC 'FUNCTI0N#PSF(0) = )i / AND THE NORMALIZED ENCIRCLED FNFRGY#') 
^«RITEa02,6065) 

6065 FORMATC'EE(INF)sI?(YESsUNO = 0)') 

REAO(101,4000)IRESP 
IFdRESP.NE.DGO TO 350 
WAVaWAVL*«DL 

DO 255 IsdNWAV 

WAVsWAV+DL 

PROD=QUEF(I)jkFF(I) 

IF(PR0D.LT.I,E-4)G0 to 255 

WAVMsWAV 

GO TO 256 

255 CONTINUE 

256 RAD*l,22*WAVM*FR 

ANGal ,2a,fcWAVM*648000,/(OP*PI) 

WRITEU02,6066)RADf ANG 

6066 FORMATC/, 'MINIMUM AIRY RADIUSa',E15,8, ' CMs%Ei5.8, 

ARC SECONDS') 

260 WRITEa02,6068) 

6068 FORMAT(/# 'ENTER MAXIMUM RADIUS TO CONSIDER(ARC SECONDS)') 
REAOdOl ,4000)RMAX 
WRITE(102,6072) 

6072 FORMAT!/, 'ENTER NUMBER OF RADIUS CALCULATIONS') 

REAOn01,4000)NR 

WRITE(102,6076) 

6076 FORMAT!/, 'ENTER DESIRED WAVELENGTH NUMBER ! 1 , 2 , , , , , 2 1 ) ' ) 
READ!lOi,4000)JW 

WRITE!102,6077) 

6077 FORMAT!/, 'CLEAR SCREEN AND RETURN TO CONTINUE') 
REAO!101,4000)OUMMY 

X=NR-i 

RCMsRMAX*DP*FR*PI/648000, 

DRaRCM/X 

DRYaRMAX/X 

DWAVaJW*! 

WAVXaWAVL+DWAV^OL 

SCAL=ENMNS*!PI*0P/!4,*FR*WAVX) )**2 
ORZaRMAX*OP<rPI / !X<kWAVX*648000, ) 

ZRa-ORZ 

YRr-DRY 

XRa-DR 


WRITE! 102, 7000 )RMAX,RCM 

7000 FORMATI'MAXIMUM RAOIUSa ' , E 1 5 . 8 , ' ARC SECONDSa ' , E 1 5 . 8 , ' CM') 
WRITE!102,7004)NR,WAVX 

7004 FORMAT! 'NUMBER OF RADIUS CALCULATIONS', 14, 's WAVELENGTHa',E15.8, 
&' CM') 

WRITE!102,7008) 

7008 FORMAT!/,' RADIUS!CM) PNT, SPREAD FN, ENCIRCLED ENERGY', 

S' !ARC SECONDS) ENERGY/SEC/CM*^^2 ' ) 

DO 340 IRrl,NR 

XRaXR+DR 

YRaYR+ORY 
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ZR=2R+0RZ 

XNUs-DNU 

DO 310 1 = 1,^100 

XNU=XNUtDNU 

X=2.i»PI<(ZR*XNU 

IF(X,GT.0,)G0 TO 285 

BES0=1. 

BESlaO. 

GO TO 300 

205 IFCX,GT,0.)G0 TO 290 
X2=X/2. 

X22=X2*X2 

X20aX22*X22 

X26*X2«*X22 

X28sX24*X20 

BESO=U-X22+X20/0.-X26/36,+X28/576,-X26*X22/10«00. 

BES0sBE80+X28*X20/51 8000, «X28*X26/2540 1600. 

X3sX22*X2 

XSsX20*X2 

X7sX26*X2 

X9=X28*X2 

BESl-X2-X3/2,+X5/12,»X7/100.+X9/28eO,-X7<fX22/86000, 
BEStsBESl+X9*X22/3628800,-X7i^X28/2032 12800, 

GO TO 300 
290 X2=X*X 
X3=X2*X 
X4 = X3i^X 
X5=X0*X 
X6sX5*X 

POsl,-0,0703125/X2+0,112152/Xa-0,5725/X6 

00=-0,125/X+0,0732a22/X3-0,227l08/X5 

8ESOe0.7978806/SQRT(X) 

AQ=X»0, 7853982 

BESOsBESO*(PO*C08(AO)-QO*SIN(AQ)) 

Pl=l ,+0.1 171 875/X2-0, 1001 956/XO+0.6765926/X6 
Q1=0,375/X-0, 1 025391 /X3+0.2775760/X5 
BES1?O,79708O6/SQRT(X) 

BQsX-2, 35619009 

BESt*BE8l+(Pl*C0S(BQ)-Ql*SIM(60)) 

300 FXO(I)aTTOT(I,JW)*BESO+XNU 
FXia)sTT0T(I,JW)+BE81 
310 CONTINUE 
A2=0. 

A3=0, 

AOaO, 

B2aO, 

B3 = 0, 

BOsO, 

DO 320 J=2, 399,2 
A2=A2+FX0(J) 

B2=82+FX1 (J) 

320 CONTINUE 
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00 330 J=3,398,4 
A3aA3+FX0(J) 

A«=A«+FXO(J-2)+FXO(J+2) 

B3=B3+FX1(J) 

B43B4-fFXl (J-2)+FXl (J+2) 

330 CONTINUE 

SAS=0NU*(l,A*A«+2.«*A3+6.-!|*A2)/a.5 

SBSaONU*(l .A^B4+2,a*83+6,4*82)/4,5 

PSF=6,*SAS/ENMN8 

8P3F=SCAL^PSF 

EE=2,*PI*ZR*SBS 

WRITEU02,7012)XR,PSF,EE,YR,SPSF 
7012 F0RMAT(5(E15.8,1X)) 

340 CONTINUE 

WRITE(1O2,6O0O) 

6080 FORMAT(/, 'DO YOU WISH TO CHANGE ANY PSF PARAMETERS? (YES= 1 ? N0=0 )' ) 
REAO(101,4000)IRESP 
IFURESP.EQ.DGO TO 260 
350 WRITE(102,6077) 

REAO(lOl,4OO0)OUMMY 

WRITE(102,6084) 

6084 FORMAT(/, 'ENTER X,Y PEL WIDTHS(ARC SECONDS)') 
READ(101,4000)XPW,YPW 
WRITE(102,6088) 

6088 FORMAT!/, 'ENTER X,Y PEL CENTER SEPARATIONSCARC SECONDS)') 
REAO(101,4000)XPO,YPO 
WRITE(I02,6092) 

6092 FORMAT!/, 'ENTER NUMBER OF PELS IN X,Y DIRECTIONS') 
REAO!101,4000)NX,NY 
DO 360 IX=1,NX 
DO 360 IY=1,NY 
360 PELB!IX,IY)sO, 

WRITE!102,6096) 

6096 FORMAT!/, 'ENTER X,Y COORDINATES OF SIGNAL STARIARC SECONDS)') 
READ!101,4000)PSX,PSY 
XNsNX 
YNsNY 

XMAXsXN*XPD 

ymax=yn*ypo 

RMAX=SQRT!!XMAX-PSX)**2'f!YMAX-PSY)**2) 

DRZ=RMAX*DP*PI/!l99,<k648000.) 

DRRxRMAX/199, 

ZRs-DRZ 

C ***** COMPUTE POINT SPREAD FUNCTION ************************** 
CALL BES!NWAV,DRZ,ZR,DNU,ENMNS,WAVL,0L) 
c ***** Integrate psf over wavelength ************************** 

IF!NWAV.GT,1)G0 to 440 
PROO=QUEF!U*FF!1)*XISL1*PI*DP*DP/4. 

DO 430 1R=1,200 
IF!KRESP.E0,1)G0 to 425 
SIGNL!IR)sPROD*PR!IR,l)/WAVL 
GO TO 430 
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425 Y=HC0K/(WAVL*8T1) 

SIGNU(lR)3pR0D*PR(IRrl)/((EXP(Y)-l.)*^AVL**6) 

430 CONTINUE 
GO TO 480 

440 MWAVcNWAV-1 
FACaXISLUPI*DP*DP*DL/4. 

DO 470 IR=1,200 
WAVsWAVL 

PROD=QUEF(l)*FF(l)*PR(IR»n 
IFCKRESP.EO.DGO TO 441 
SUMsPROD/WAV 
GO TO 442 

441 Y=HCOK/(WAV*STJ) 

SUM=PROD/(CEXP(Y)»l ,)!»WAV**6) 

442 IF(NWAV,EQ,2)G0 TO 460 
DO 450 JW32»MWAV 
WAVsWAV+DL 

PROD=QUEF(JW)*FF(JW) 

IF(PR0D.LT,I,E-4)G0 to 450 
IF(KRESP.EQ.1)G0 TO 446 
SUMaSUM-fPROD*PR(IR#JW)/WAV 
GO TO 450 

446 Y=HCOK/(WAV*STl) 

SUM=SUM+PR00*PR(IR,JW)/((EXP(Y)-1.)aWAV**6) 

450 CONTINUE 
460 ^AVaWAVtDL 

PROOsQUEF(NWAV)*FF(NWAV)*PR(IRrNy^AV) 

IFCKRESP.EO.DGO TO 465 
SUMsSUM+PROD/WAV 
60 TO 466 

465 Y=HCOK/(WAV*STD 
8UMsSUM+PR00/((EXP(Y)-l.)*NAV**6) 

466 SI6NUIR) = SUM*FAC 
470 CONTINUE 

C ***** INTEGRATE SIGNAL OVER PELS ****************************** 
480 DPXsXPW/20, 

DPYsYPW/20, 

DO 500 IXsl#NX 
DO 500 IYs 1,NY 
PEL3(IX,IY)aO. 

XlsIX-l 

YIsIY-1 

XXs-DPX/2.+XI*XPD 

SUM=0, 

DO 490 IDXSD20 

YYs*DPY/2.+YI*YPD 

XXsXXtDPX 

DO 490 IDY=D20 

YYaYYtDPY 

DIF=(PSX-XX)<k*2+(PSY-YY)**2 

DST=SQRT(0IF) 

FRACsDST/DRR 
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NUMsFRAC 

IF(NUM,GE.199)G0 TO 500 
XNUM=NUM 
FRAC=FRAC-XNUM 
F A C® 1 « 

IF((IDY.EQ,U.OR,(IDY.EO,20))FAC = 0,5 
IF((IDX.EQ,l),OR.(IDX.EQ.20))FAC=FAC*0,5 

SUMaSUMfFAC*(SIGNL(NUM+l)+FRAC*(SIGNL(NUM+2)-SIGNL(NUM+n ) ) 

490 CONTINUE 

PELS ( I X,lY)aSUM*ENMNS*PI*DPX*DPY* (PI *DP/ 1296000. )**2 
500 CONTINUE 

C ***** CALCULATE DIFFRACTION EFFICIENCY *********************** 

WAV=WAVL"OL 
DO 920 NWal,NWAV 
WAVsWAV+OL 

PRODsQUEF(NW)*FF(NW) 

IF(PR0D.LT,1,E-4)G0 TO 920 

DO 920 IX=1,NX 

DO 920 IY=1,NY 

XI=IX-1 

YisIY-i 

XX=-0PX/2,+XI*XPD 

SUMsO. 

DIFE(IX,IY)»0, 

DO 900 IDX=l,20 

YY=«DPY/2.+YI*YP0 

XXsXX+DPX 

DO 900 IDY=1,20 

YYsYY+DPY 

DSTsSQRT((PSX-XX)**2+(PSY-YY)**2) 

FRACsDST/ORR 

NUMsFRAC 

IF(NUM.GE,199)G0 TO 9J0 
XNUMsNUM 
FRACsFRAC-XNUM 
FAC=1. 

IFCaDY.EQ.n,0R,(IDY.EQ,20))FACs0,5 

IF((IDX,EQ,n,OR,(IDX.EQ,20))FACsFAC*0,5 

8UMsSUM+FAC*(PR(NUM+l ,NW)+FRAC*(PR(NUM+2,NW)-PR(NUM+1 ,NW) ) ) 

900 CONTINUE 

DIFEdX, IY)=SUM*ENMNS*PI*DPX*DPY*(PI*DP/(WAV*1296000. ) )**2 
910 V«RITE(l0a,70t7)WAVdX,IY»DIFE(IX,IY) 

7017 F0RMAT('WAVs%E 15.8, ' IXs'dlr' IY=%Ilf' OIFEs' , E 15 ,8) 

920 CONTINUE 

C CALCULATE BACKGROUND SIGNAL **************************** 

WRITE(102,6100) 

6100 FORMAT(/, 'ENTER COSMIC BACKGROUND VISUAL MAGNITUDE') 
REAO(101,4000)XMB 
X1MBsO.37353*1O.O**(-O.4*XM0-5.) 

IFCKRESP.EQ.UGO TO 510 
XIMLBs(XIM8*10,0**5)/(l,0682*HC) 

GO TO 515 
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510 WRITE(102,6l0a) 

610^ FORMATC/, 'ENTER BACKGROUND EFFECTIVE TEMPERATURE (KEUVIN)' ) 
READ(101,4000)BEFT 
CAU PHOI(BEFT»HCOK,DLA,HC,XINT) 

XIMLB?XIMB*XINT 

C ***** INTEGRATE OVER WAVELENGTH *********************:**★****** 
515 IFCNWAV.GT.DGO TO 520 
IF(KRESP.EQ.1)G0 TO 517 
SUM=XIMLB*OUEF(l)*FF(l)*WAVL 
GO TO 550 

517 Y=HCOK/(WAVL*BEFT) 

SUMsXIMU8*QUEF(l)*FF(l)/((EXP(Y)-l,)*WAVL**a) 

GO TO 550 

520 WAVsWAVL 
IF(KRESP.EQ.1)G0 to 521 
SUMaQUEF(l)*FF(l)*WAVL 
GO TO 522 

521 Y=HC0K/(WAV*8EFT) 

SUMcQUEF(1)*FF(1)/((EXP(Y)-1,)*WAV**0) 

522 IF(NWAV.EQ,2)G0 TO 5^0 

mwav=nwav-i 

DO 530 I=2,MWAV 

wav=wav+ol 

PRODsQUEF(I)*FF(I) 

IF(PR00.LT.I,E-0)G0 TO 530 
IF(KRESP.E0.1)G0 to 528 
SUM=SUMfPROO*WAV 
GO TO 530 

528 Y=HC0K/(WAV*BEFT) 

SUMsSUM+PR0D/((EXPCY)-1,)*WAV**4) 

530 CONTINUE 
5^0 WAVaWAV+DL 

IFCKRESP.EQ.DGO to 5«5 
SUMaSUM+QUEF(NWAV)*FF(NWAV)*WAV 
GO TO 546 

545 YsHC0K/(WAV*BEFT) 

8UMaSUM+QUEF(NWAV)*FF(NWAV)/( (EXP(Y)-l ,)*WAV**4) 

546 SUMaSUM*OL*XIMUB 

550 BGRN0a0,7854*DP*DP*SUM*XPW*YPW 
IRESPaO 

WRITE(102,6106) 

6108 FORMAT(/,'D0 YOU WISH TO INCLUDE A BACKGROUND STAR? ( YESal ; NOaO ) ' ) 
READ(101,4000)IRESP 
IFCIRESP.NE.DGO TO 700 
WRITE(102,6112) 

6112 FORMATC/, 'ENTER VISUAL MAGNITUDE OF BACKGROUND STAR') 
REA0(101,4000)XVN 
XIMV=0,37353*10.0**(-0,4*XVM-5.) 

IF(KRESP.EG.1)G0 TO 560 
XIMLVa(XIMV*10,0**5)/(l,0682*HC) 

GO TO 570 

560 WRITE(102,6116) 
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6U6 FORMAT(/, 'ENTER BACKGROUND STAR EFFECTIVE TEMPERATURE (KELVIN) ' ) 
READ(101,4000)ST2 
CALL PH0I(ST2,HC0KfDLA,HC,XINT) 

XIMLV=XIMV*XINT 

570 WRITE(i02,6120) 

6120 FORMAT!/, 'ENTER X,Y COORDINATES OF BACKGROUND STAR(ARC SECONDS)') 
REA0a0i,/l000)BSX,8SY 
DXs8SX-XMAX 
DY=8SY-YMAX 
D00a8SX*8SX+BSY*8SY 
0X0s0X*DX+8SY*eSY 
DXY=DX*DX<fDY*DY 
OOYsBSX*8SX+DY*OY 
OMAXpO, 

DMINslOOOO, 

IF(000,GT,DMAX)OMAXsOOO 

IF(OXO,GT.DMAX)OMAXsDXO 

IF(OXY,GT,OMAX)DMAXrOXY 

IF(D0Y,GT,0MAX)DMAX=D0Y 

IF(000,LT,OMIN)OMINsDOO 

IF(DX0.LT.DMIN)DMIN=DX0 

IF(DXY,LT,0MIN)DMIN=DXY 

IF(OOY,LT.OMIN)DMINcDOY 

OMAXpSQRT(OMAX) 

OMINaSQRT(DMIN) 

IF(0X,GT,0,)GO TO 580 
IF(DY,GT.O.)GO TO 600 
DMINpO, 

GO TO 610 

580 IF(DY.GT.O.)GO TO 610 
DMINsOX 
GO TO 610 

600 IF(0X,6T,0.)GO TO 610 
DMINbOY 

610 RMAXsOMAX-DMIN 

DRZ=RMAX*DP*PI/(199,*6«8000.) 

DRRaRMAX/199, 

ZRsDMIN*DP*PI/6a8000.-DRZ 

CALL BES(NWAV,DRZ,ZR,DNU,EMMNS,WAVLrDL) 

***** INTEGRATE OVER WAVELENGTH ****************************** 

IFCNWAV.GT.DGO to 630 

PROObQUEF ( 1 ) *FF ( 1 ) *ximlv*pi*dp*dp/^, 

DO 620 IRsl,200 
IFCKRESP.EQ.DGO to 615 
SIGNL(IR)=PR0D*PR(IR,1)/WAVL 
GO TO 620 

615 Y=HC0K/(WAVL*ST2) 

SIGNL(IR)=PR00*PR(IR,1)/((EXP(Y)-1.)*WAVL**6) 

620 CONTINUE 
GO TO 670 

630 MWAV=NWAV-1 

FACsXIMLV*PI*DP*OP*DL/«, 



DO 660 IRslf200 
WAVsWAVU 

PRODsQUEF(l)*FF(n*PR(IRf 1) 

IFCKRESP.EQ.l)GO tO 631 

SUMaPROO/WAV 

GO TO 632 

631 YaHC0K/(WAV*ST2) 

SUM8PR0D/((EXP(Y)-1.)*WAV**6) 

632 IF(NWAV.EQ,2)GO TO 650 
DO 6A0 Jrta2,MWAV 
WAVaWAVfDL 

PRODaQUEF(JW)^FF(JW) 

IF(PR0D,LT.l,E-4)60 TO 6^0 
IFCKRESP.EQ.DGO TO 636 
SUMaSUM+PROD*PR(IR,JW)/WAV 
GO TO 6A0 

636 Y=HC0K/(WAV*ST2) 

SUMaSUM+PROD^PRdR, JW)/((EXP(Y)-1.)*WAV*^6) 

6a0 CONTINUE 
650 WAVaWAV+DL 

PROD=QUEF(NWAV)*FF(NWAV)*PR(IR,NWAV) 

IFCKRESP,EQ.UGO TO 655 
SUMaSUM+PROO/WAV 
GO TO 656 

655 Y = HCOK/(WAVJkST2) 

SUMaSUM+PR0D/((EXP(Y)*l,)^WAV**6) 

656 SIGNU(IR) = SUM<rFAC 
660 CONTINUE 

c ***** Integrate signal over pel ******************************* 

670 DPXaXPW/20. 

DPYaYPW/20, 

DO 690 IX=1#NX 
DO 690 lYalfNY 
PELB(IX,IY)=0. 

XIaIX-1 

YX=IY-1 

XXa-0PX/2,+XI*XP0 

SUMaO. 

DO 680 IDXal,20 
YYa-DPY/2. + YI<(YPD 
XXaXX+DPX 
DO 680 IDY=1#20 
YYaYY+DPY 

DIF=(BSX-XX)Ai*2t(BSY-YY)**2 

DSTaSQRT(OIF)-DMIN 

FRACbDST/DRR 

NUMrFRAC 

IF((NUM,GE.199).OR.(NUF.LT.O))GO TO 690 
XNUMaNUM 

fracsfrac-xnum 

FACsl, 

lFC(IDY,EQ.l),0R.(I0Y.EQ,20))FACa0,5 
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IF((IDX,EQ,l),OR, (IDX.EQ.20))FAC=FAC*0,5 

SUM=SUM+FAC*SIGNL(NUM+l)+FAC*FRAC*(3IGNL(NUM+2)-SIGNL(NUM+n ) 
680 CONTINUE 

PELS (IX, I Y)=SUM*ENMNS*PIaDPX*OPY*CPI*DP<^FR/ 1 296000, )**2 
690 CONTINUE 

C CALCULATE DIFFRACTION EFFICIENCY ************^;******a*** 

WAV=WAVL-DL 
DO 950 NWs1,nV^AV 
WAVsWAVtOL 

PROOsQUEF(NW)*FF(NW) 

IF(PROD.LT,l ,E-a)GO TO 950 

00 950 IX=1,NX 

00 950 IY=1,NY 

XI=IX-I 

YI=IY-1 

XXs-0PX/2,+XH^XPD 

SUMsO, 

0FE(IXrIY)=0, 

DO 930 IDX=l,20 

YY=-DPY/2.+YI*YPD 

XX=XX+DPX 

DO 930 IDY=1,20 

YY=YY+DPY 

DST=SQRT((BSX-XX)**2f (BSY-YY)**2)-DMIN 

FRAC=DST/ORR 

NUMsFRAC 

IF((NUM.GE,199).0R,(NUM,LT,0))G0 to 940 

XNUMsNUM 

FRACsFRAC-XNUM 

FAC®1 

IF((IOY,EQ,l),OR.(IDY.EQ,20))FAC=0,5 

IF((I0X.EQ,l).OR.(IDX.EO.20))FAC=FAC*0,5 

SUM=SUM+FAC*(PR(NUM+l,NW)+FRAC*(PR(NUM+2,NW)-PR(NUM+!,NW))) 

930 CONTINUE 

DFECIXf lY)s SUM*ENMN3ikPI*DPX^DPY* (PI *DP/(WAV*1 296000. ) )**2 
940 WRITE(102,7027)WAVf IX,IY,0FE(IX,IY) 

7027 F0RMAt('WAVa',E15,8,' IXaMl,' IY=M1,' DlFE=',E15,8) 

950 CONTINUE 

700 WRITE(102,6124) 

6124 FORMAT!/, 'ENTER RMS READOUT NOISE(ELECTRONS/PEL) ') 
READ(101,4000)RON 
RN=RON*RON 
!ARITEa02,6128) 

6128 FORMAT!/, 'ENTER DETECTOR TEMPER ATURE ( KELV IN )' ) 
READ(101,4000)OTEM 
IF(DTEM,LT,4,)G0 TO 710 

FAC=*6705,487/DTEM+4.0737*DTEM/(U08,+DTEM) 

FAC-EXP(FAC)*(29,3E+9) 

DC=FAC'kDTEM**l,5 
GO TO 720 

710 WRITE(i02,6132) 

6132 FORMAT!/, 'ENTER MEAN DARK CURRENT (ELECTRONS/PEL/SEC) ' ) 
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READ(101,4000)DC 

720 WRITEO02,6136) 

6136 FORMAT(/, 'ENTER OBSERVATION START AMO END TIMES ( SECONDS) ' ) 
READ(lOl,aoOO)STRT#ENDT 
WRITEdoa, 61410) 

61410 FORMAT(/, 'ENTER NUMBERC>1) OF TIME CALCULATIONS') 

READ(10l,4|000)NT 

X=NT-l 

DT=(ENDT-STRT)/X 
WRITE(102, 614140 

614IA FORMATC/, 'ENTER START AND END SIGNAL TO NOISE RAT lOS ( S/N) ' ) 
READ(l0l,4|000)SNSTfSNND 
WRITE(102, 61418) 

61418 FORMAT!/# 'ENTER NUMBEROl) OF (S/N) CALCULATIONS') 
REAO(10l#4|000)NS 
'Y*NS*1 

DNS=(SNNO-SNST)/Y 

WRITEU02,6077) 

REAO(101,4»oOO)OUMMY 

725 WRITE(102,61S2)DP 

6152 FORMAT(/#'TELESCOPE ENTRANCE PUPIL DIAMETER ',F7.3#' CM') 
WRITE(l02,6156)FR 

6156 FORMAT('OPTICAL SYSTEM FOCAL RATIO '#F5,2) 

WRITE(102,6160) 

6160 F0RMAT(/# 'SOURCE STAR CHARACTERISTICS') 

WRITE(102,6164OXSl 

61641 FORMATCSX, 'SOURCE MAGNITUDE '#F5,2) 

IF(KRESP,EQ,0)60 TO 726 
WRITE(102,6168)ST1 

6168 FORMATCSX, 'SOURCE TEMPERATURE(KELVIN) ',F9.3) 

726 WRnEU02#6l72) 

6172 FORMAT(/# 'COSMIC BACKGROUND CHARACTERISTICS') 

WRITE(102#6176)XMB 

6176 FORMATCSX# 'COSMIC BACKGROUND MAGNITUDE '#F5.2) 

IF(KRESP,EO.O)GO TO 727 
WRnE(l02#6ie0)BEFT 

6180 F0RMAT(5X, 'COSMIC BACKGROUND TEMPERATURE(KELVIN) '#F9.3) 

727 IFURESP.NE.DGO TO 730 
WRITE(l02#6184O 

61841 FORMAT!/# 'BACKGROUND STAR CHARACTERISTICS') 

WRITE(102#6188)XVM 

6188 F0RMAT(5X, 'BACKGROUND STAR MAGNITUDE '#F5.2) 

WRITE!102#6192)ST2 

6192 F0RMAT!5X# 'BACKGROUND STAR TEMPERATURE(KELVIN) '#F9.3) 

730 WRITE(102#6196) 

6196 FORMAT!/# 'SYSTEM SPECTRAL CHARACTERISTICS') 

WRITE(102#6200) 

6200 FORMAT (4»X# 'WAVELENGTH (CM) '#SX#'QUAN, EFFIC.',3X# 'FILTER FUNCTION') 
WAVsWAVL 
DO 7A0 Ial#NWAV 
PROOsQUEF(I)*FF(I) 

IF(PR0D,LT.1.E«A)60 TO 739 
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WRITE (102,6304) WAV #QUEF (I), FF( I) 

6204 F0RMAT(3(6X,E11,4)) 

739 WAVaWAV+DL 

740 CONTINUE 
WRITE(102,6208) 

6208 FORMAT(/, 'DETECTOR CHARACTERISTICS') 

WRITE(1Q2,6212)XPW,YPW 

6212 F0RMATC5X, 'PEL DIMENSIONS (ARC SECONDS) :Xs*,F 8. 4, ' Ys',F8,4) 

WRITE(102,6216)XPD,YPD 

6216 F0RMAT(5X, 'PEL CENTER SEPARATION(ARC SECONDS) : X= F8, 4 , ' Y=', 

&F8,4) 

WRITE(102,6220)PSX,PSY 

6220 F0RMAT(5X, 'SOURCE STAR COORDINATES(ARC SECONDS) : Xs F8.4 , ' Y=', 

&F8,4) 

IFdRESP.NE.DGO TO 750 
WRITE(1Q2,6224)8SX,BSY 

6224 FORMAT(5X, 'BACKGROUND STAR COORDINATES (ARC SECONDS )! X= ', F8 , 4 , 

&' Ya'iF8.4) 

750 WRITE(102,6228)NX,NY 
6228 F0RMAT(5X, 'ARRAY SIZE:Xs',I2, ' BY Y=',I2) 

WRITE(102,6232)R0N 

6232 F0RMAT(5X,'RMS READOUT NOISE(ELECTRONS/PEL) : ',F5,2) 
IF(DTEM.LT,4.)G0 to 760 
WRITE(102,6236)DTEM 

6236 F0RMAT(5X, 'DETECTOR TEMPERATURE(KELVIN) : ',F8.4) 

760 WRITE(102,6240)DC 

6240 F0RMATC5X, 'DARK CURRENT(ELECTR0N8/PEL/SEC0N0) : ',£15,8) 
WRITE(102,6246) 

6246 FORMAT(/, 'SIGNAL TO NOISE RATIO AS A FUNCTION OF OBSERVATION', 

&' TIME*) 

WRITE(102,6250) 

6250 FORMAT('AND PEL NUMBER') 

WRITE(102,6254) 

6254 F0RMAT(4X,'TIME(SEC0NDS) ',2X,' (S/N) ',2X,'PEL NUMBER') 

DO 770 IXal^NX 
DO 770 IYal,NY 
T=STRT-OT 
DO 770 ITx1,nT 
T=T+DT 

XNUMspELS(IX,IY)*T 

X0N0M=XNUM+8GRND*T+PELB(IX,IY)*TtRN+DC*T 

SNpsXNUM/SQRT(XDNOM) 

WRITE(102,6256)T,SNP,IX,IY 
6256 F0RMAT(5X,E12,5,6X,E11.4,2X,213) 

770 CONTINUE 

WRITE(102,6260) 

6260 F0RMAT(/,5X,'(S/N)',4X,' SECONDS ',2X,'PEL NUMBER') 

DO 780 IXal,NX 
DO 780 IY=1,NY 
SNsSNST-DNS 
DO 780 13=1, NS 
SNsSNfONS 
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A=SN*SN/(PELS(IX,IY)*a2) 

B=A*(PELS(IX,IY)+BGRND+PELBaXrIY)+DC)*0.5 

C3A*R0N 

n=B+SQR7CBAB+C) 

WRITE(l02,626a)SN,Tl,IX,IY 

626^ F0RMAT(5X,F5,2,5X,E12.5,2X,2I3) 

780 CONTINUE 

WRITE(lQ2,626e) 

6268 FORNAT(/»'DO YOU WISH TO CHANGE SOURCE MAGN I TUDE? ( YES= 1 ; N0=0 ) ' ) 
READ(l01,aOOO)JRESP 
IF(JRESP.EQ.O)GO TO 800 
WRITEa02f6272) 

6272 FORMAT!/, 'ENTER NEW SOURCE MAGNITUDE') 

READ(101#^000)XSM 

XNEWs0.37353*10,0**(-0,4*XSM) 

XNEWsXNEW/XIMS 
DO 790 IXsl,NX 
DO 790 IYsl,NY 

790 PELSUX,IY)=PELSCIX,IY)AXNEW 
800 WRITE(102,6276) 

6276 FORMAT!/, 'DO YOU WISH TO CHANGE COSMIC BACKGROUND MAGNITUDE?', 
&'!YESs1?N050)') 

READ!lOl,AOOO)JRESP 
IF!JRESP.EQ.O)GO TO 810 
WRITE!102,6280) 

6280 FORMAT!/, 'ENTER NEW COSMIC BACKGROUND MAGNITUDE') 
READ!l0l,«OOO)XMB 
XNEW*0.37353<^10.0*^!-0,«*XMB) 

XNEW#XNEW/XIMB 
BGRNDsBGRND*XNEW 
810 WRITE!102,628A) 

6289 FORMAT!/, 'DO YOU WISH TO CHANGE BACKGROUND STAR MAGNITUDE?', 
&'!YES=1 ;NO=0) ') 

READa01,9000)JRESP 
IF!JRESP.EQ,0)G0 TO 830 
WRITE!102,6288) 

6208 FORMAT!/,'ENTER NEW BACKGROUND STAR MAGNITUDE') 
READ!101,9000)XMV 
XNEW=0,37353^10,0<^*!-0.9*XMV) 

XNEW5XNEW/XIMV 
DO 820 IX=1,NX 
DO 820 IY»1,NY 

820 PELBaX,IY)=PElB!IX,IY)*XNEW 
830 WRITE!102,6292) 

6292 FORMAT!/, 'DO YOU WISH TO CHANGE DETECTOR EFFECTS?!YES=1 ,NOsO) ') 
READaOl,9000)JRESP 
IFURESP.EQ.OGO TO 700 
WRITEC102,6296) 

6296 FORMAT!/, 'DID YOU MAKE ANY CHANGES? !YES=1 jNOsO) ') 
READU0U9000)JRESP 
IFIJRESP.EQ.DGO TO 725 
STOP 
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END 

SUBROUTINE PHO I C TENP, HCOK , DLA , HC , X INT ) 
implicit double PRECISI0N(A-H,0-Z) 

COMMON /PH0F/XK(i»0) 

X=3.8E*5 

Y=HC0K/(X*TEMP) 

SUMsXKCl )/((EXP(Y)-l,)*X**5) 

x=x+dla 

DO 10 IK=2,38 

ICHEK=2*(IK/2) 

FAC52, 

IF(IK.EO.ICHEK)FAC=4. 

YsHCOK/(X*TEMP) 

SUMsSUM+FAC*XK(IK)/((EXP(Y)-l,)*X**5) 

XsX+DLA 
10 CONTINUE 

YsHCOK/(X*TEMP) 

SUMsSUM+XK ( 39) /( (EXP (Y)-1.)*X**5) 

XINT=3./(SUM*DLA*HC) 

RETURN 

SUBROUTINE BES ( NWAV r DRZ, ZR, DMU, ENMNS f WA VL r DL ) 

IMPLICIT DOUBLE PRECISI0N(A-H,0-Z) 

COMMON /aESF/FXO(^00),PR(^OOf2i),TTOT(^00,^l),QUEF(^l),FF(2l) 

TW0PI=6. 283185307 

WAVsWAVL-OL 

DO 60 NW=l,NWAV 

WAVaWAV+DL 

OR=DRZ/WAV 

Z=ZR/WAV 

PROO?QUEF(NW)*FFCNW) 

DO 60 IR=1,200 
PR(IR,NW)=0. 

Z = Z-^DR 

IF(PR00.GT.1.E-4)G0 TO 5 
GO TO 60 
5 XNUs-DNU 

DO 31 1=1,400 
XNU=XNU+DNU 
X=TW0PI*Z*XNU 
IF(X,GT.O,)GO TO 10 
BS=1. 

GO TO 30 

10 IF(X,GT.4,)G0 to 20 
X2=X/2, 

X22=X2*X2 

X24=X22*X22 

X26=X24*X22 

X26=X24*X24 

BS=1 ,-X22+X24/4,-X26/36,+X28/576.-X28*X22/ 14400, 
8S=BS+X28*X24/51 8400. -X2fi*X26/2540 1600, 

GO TO 30 
20 X2=X^X 
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X3sX2*X 

X4sX3*X 

X5sX4*X 

X65XS*X 

Psl, *0,0703 125/X2+O.U2152/X4-0,5725/X6 

Qs*0,125/X+0,0732422/X3-0,227108/X5 

BSaO,7978846/SORT(X) 

AQaX-0. 7853982 
8S=8S*(P*C0S(AQ)-Q*SIN(AQ)) 

30 FXOiI)=TTOT(I,NW)*BS*XNU 

31 CONTINUE 
A2=0, 

A3“0 • 

A4a0, 

00 40 Ia2,399,2 
A25A2+FX0(I) 

40 CONTINUE 

00 50 Ia3,398,4 
A3aA3+FX0(I) 

A4aA4+FX0(I*2)tFX0(I+2) 

50 CONTINUE 

8ASa0NU*(1.4*A4f2,4*A3+6.4*A2)/4.5 
PR(IR,NW)ae,*SAS/ENMNS 
60 CONTINUE 
RETURN 
END 
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